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Charge particle transport across nanogaps is studied theoretically within the Schrodinger-Poisson mean field 
framework and the existence of limiting current investigated. It is shown that the choice of a first order 
WKB wavefunction as the transmitted wave leads to self consistent boundary conditions and gives results 
that are significantly different in the non-classical regime from those obtained using a plane transmitted wave. 
At zero injection energies, the quantum limiting current density ( J c ) is found to obey the local scaling law 
J c ~ V g / D 5 ~ 2a with the gap separation D and voltage V g . The exponent a > 1.1 with a — > 3/2 in the 
classical regime of small de Broglie wavelengths. These results are consistent with recent experiments using 
nanogaps most of which are found to be in a parameter regime where classical space charge limited scaling 
holds away from the emission dominated regime. 



The current-voltage characteristics across nanogaps 
is a subject of much interest across a wide variety of 
fields. An important quantity is the maximum current 
that can be transmitted across the gap. In classical 
physics, the limit exists due to the mutual repulsion be- 
tween the charged particles and the current is said to 
be space charge limited. In a mean field picture, the 
mutual repulsion gives rise to a potential barrier. As 
the current reaches the limiting value, the barrier height 
increases and the transmitted particles barely cross the 
barrier. Beyond this value, the barrier height oscillates, 
some electrons get reflected and steady state transport 
is no longer possible. The Child-Langmuiri^ law and its 
many generalizations^— serve well to predict the max- 
imum transmitted current in macro devices where the 
available phase space volume is large and quantum ef- 
fects can be safely ignored. 

In nanostructures however, the small physical size to- 
gether with low applied voltages or injection energies 
pushes the average per electron phase space volume to- 
wards the plank cell limit. Thus, in addition to space 
charge, quantum effects must be considered in exploring 
the existence of such a limit in nanoelectronics, vacuum 
microelectronics or devices such as the scanning tunnel- 
ing microscopes. In particular, the possibility of tunnel- 
ing leads to the question: is there a limit on the maximum 
transmitted current density in the quantum mechanical 
case ? Theoretical studies^ suggest that quantum tun- 
neling pushes up this limit and by orders of magnitude 
in the low injection energy regime where quantum ef- 
fects should be considerably more important. In fact, 
within the framework considered in Ref. [(| the limit can 
indeed be very large in the very low energy regime as 
we shall show here. This singular behavior, signalling a 
sharp departure from the classical prediction, warrants a 
fresh look at the basic assumptions involved, especially 
since, many of them are retained in more sophisticated 
theories^— which take into account the fermionic nature 
of electrons by incorporating the exchange-correlation 
potential. 



In its simplest form, the mean field quantum (Hartree) 
formalism involves solving the coupled Schrodinger and 
Poisson equations 



h 2 cPtp 



2m dx 2 
d 2 V _ eH 2 



eV(x)ip = Evb 



dx 2 



(1) 

(2) 



in the region [0, D] where D is the size of the gap, e is the 
magnitude of the electronic charge, E is the energy of the 
electron and V(0) = while V(D) = V g . Such a model 
is expected to hold when the electron density in the gap 
is sufficiently high that their mutual interaction must be 
accounted for, but is low enough to neglect the effects of 
the exclusion principle. It is assumed here that V is time- 
independent and there is a steady current flowing across 
the gap with a current density J — eih(ip*ip —ipip* )/2m. 
The boundary conditions or initial values for ip are to be 
determined takinginto account this fact. 

Following Ref. |6j, we write the wavefunction ip as 
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where r(x) and 0(x) are real, E = E/eV s is a dimension- 
less energy and n s is a characteristic density. In terms 
of other dimensionless quantities x — x/D, V — eV/E, 
J = J/J s , n = n/n s = \ijj\ 2 /n s and 4> g = eV g /E, 
where V s = h 2 /(2meD 2 ), J s = e h 3 /{4:m 2 eD 5 ), n s = 
eoh 2 1 (2me 2 D A ) , the coupled Schrodinger and Poisson 
equations can be expressed as 
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where A = 2J/E 3 / 2 is a dimensionless perveance. Note 
that once r(x) is known, 9(x) can be determined indepen- 
dently with an arbitrary phase (9(1) = 0) at the bound- 
ary x = D. 

Assuming now that at x ~ D, t)j(x) = Ce lv{ - D ^ x l h , it 
follows that r(l) = (A/4)V 2 /(l + ff )i/4 while r'(l) = 0. 
Here p(a;) = y/2m(E + eV(x)) is the classical momen- 
tum. Note that the plane wave assumption above pre- 
supposes a constant potential immediately beyond the 
domain of interest^ irrespective of the nature of the in- 
terface at x = D, the length scales involved and the self 
consistent nature of the problem. 

At a given scaled injection energy E, the dimensionless 
perveance A is increased till no solution exists. This gives 
the critical current density J c = J S A C E 3 / 2 /2. In Ref. |g, it 
has been observed that for small values of E, A c exceeds 
the classical values by orders of magnitude. 

The behaviour of the system at low injection energies 
(E) can be analysed by neglecting the second term alto- 
gether in Eq. (Q} at finite values of A. It is then easy to 
see that at very low energies, a solution exists as A is in- 
creased atleast so long as the second term in Eq. (0} can 
be neglected. This behaviour is however peculiar to the 
initial conditions specified above. We show later that the 
same equations albeit with a different initial condition for 
r'(l) imposes a limit on A beyond which no solution exists 
even at low injection energies. 

Apart from the E — > behaviour, the semiclassical 
consistency of the transmitted plane wave assumption is 
also worth investigating. If the potential V for a particu- 
lar value of E and A is well behaved in the region [0, D], it 
is reasonable to expect that away from classical turning 
points of the potential V, the transmitted wavefunction 
should have the standard first order WKB form: 



iP sc (x) = -^=eif*^ dx ' +i v (7) 
Vp( x ) 

as x approaches D. Here p(x) is the classical momentum 
while C is assumed to be real with all the phase infor- 
mation dumped in if. Assuming this form to be true at 
x = xq < D, it is possible to compare r(x) and 9(x) 
with the semiclassical amplitude and phase by integrat- 
ing backwards from xo = xq/D. Thus, starting with 
r(x ) = r sc (x ) = C/y/p(x ) and 9(x ) = 9 sc (x ) = 
J p(x')dx' /H+tp, one can compare r(x) and 6(x) for val- 
ues of x < xq with the semiclassical (WKB) predictions: 
r sc (x) = Cj \Jp(x) and 9 sc (x) = f x p(x')dx'/h + ip = 
9(xq) — p(x')dx' jh. A comparison of the amplitudes 
for E = 1 and 1000 is shown in Fig. 1 for A = 15 and 
4> g = 0.5. While 9(x) and 9 sc (x) agree reasonably at 
E = 1 near xq — 0.9 (not shown here), the amplitude r(x) 
has very little agreement with r sc (x) even slightly away 
from x = Xq (see Fig. 1). At E = 1000 however, the am- 
plitude oscillates about the the semiclassical prediction 
while 9(x) and 9 sc (x) are practically indistinguishable in 
the entire region. 
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FIG. 1. (Color online) A comparison of the amplitude r(x) 
obtained using the plane wave boundary condition with the 
semiclassical prediction r sc (x) (dashed curves). The agree- 
ment is poor at E = 1 but improves for E — 1000. Here 
A = 15 and <f> g = 0.5 while xq — 0.9. 

These results are indeed not very surprising since the 
plane wave approximation is only the first term in a semi- 
classical expansion of tp(x) — e^ s ^: S(x) — Sq(x) + 
iS^x) + (i) 2 S 2 (x) + ... where S (x) = f x p(x')dx>. By 
dropping Si and subsequent terms in a plane wave expan- 
sion, it is only to be expected that the amplitude is not 
reproduced accurately even for large E while the phase 
information is well approximated as E increases. 

For a well behaved potential, a better approxima- 
tion for the transmitted wave near x — D should be 
the first order WKB wavefunction in Eq. §7§, while for 
x ~ 0, it should be a superposition of a right and a 
left moving wave of the WKB form: ip S c(x — 0) = 
(A/^/pJ^))eif a, P^"> dx ' + (B/^/pJx))e-a !C ^ x ^ dx ' . 
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FIG. 2. (Color online) A comparison of the amplitude r(x) 
obtained using the new initial condition with the WKB pre- 
diction r ac (x) for A = 15 and <j> g — 0.5. The agreement is 
better even at E = 1. At E = 100, there is practically no 
reflection. The agreement is therefore good even near x = 0. 
The solid lines are the exact results; the dashed/dotted lines 
are the WKB predictions. 

Assuming the semiclassical form (Eq. (J7J|) to be valid 
at x = D, the initial conditions for r(x) are r(l) = 
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(A/4) 1 /2/(l + (fig) 1 / 4 and r'(l) = -\/'(l)A 1 / 2 /[8(l + 
</> 9 ) 5/4 ]. Note that r'(l) depends on This does 

not preclude r'(l) = 0; rather it forces the potential to 
assume a constant value smoothly at x ~ D in order to 
be semiclassically consisten t 10 ' 11 . In practice, Eqns. (U]) 
and (O are solved as an initial value problem starting at 
x = 1 and integrating backwards till x = by choosing a 
value for such that V(0) = 0^. For the parameter 

values studied by us however, V'(l) ^ for any allowed 
solution set. Not surprisingly, a comparison of r(x) and 
r sc {x) shows a much better agreement now as shown in 
Fig. 2. 

Imposition of the semiclassical boundary condition 
also removes the singular behaviour at low scaled energies 
(E — > 0). Neglecting the second term in Eq. (Q}), and 
using the new initial conditions, the amplitude equation 
gives r(x) — ax+b where a+b = (X/4) 1 / 2 / (l + cfig) 1 / 4 and 
a = -l/'(l)A 1 / 2 /[8(l+</> 9 ) 5/4 ]- Inserting r(x) in the Pois- 
son equation (Eq. [5]) and using the boundary condition 
V(0) = gives V(x) = a 2 x 4 /12 + abx 3 /3 + b 2 x 2 /2 + ex. 
Finally on using V(l) = (f> g and demanding that r'(l) = 
-(a 3 /3 + ab + b 2 + c)A 1 / 2 /[8(l + (figf'% a real solution is 
found to exist only when (48£ 6 ) 2 - (96£ 6 + 12£ 2 (j> g )\ - 
A 2 /2 > where E = (1 + (fig) 1 / 4 . This imposes an upper 
limit on the perveance A at low injection energies which 
matches with the critical perveance A c at small E ob- 
tained by solving the full set of equations (Eqns. (j4j and 
(0)) for different values of A. For instance, the inequal- 
ity above predicts that a real solution exists at <p g = 0.5 
for A < 38.33 at very small values of E. This is a good 
approximation to the actual value at small E as can be 
seen in Fig. 3 where A c is plotted for different scaled in- 
jection energies for <fi g — 0.5, and -0.5. Note that A c 
does not increase by orders of magnitude even at low in- 
jection energies. Rather the quantum regime manifests 
itself differently at various applied potentials. While A c 
is about double the classical value at <fi g = 0.5, it is less 
than the classical value at <fi g = —0.5 where reflection 
(rather than tunneling) dominates. In general, at most 
injection energies, the maximum current evaluated quan- 
tum mechanically falls short of the classical prediction. 

When the injection energy is exactly zero, (E = 
in Eq. (fTj)). the semiclassical formalism in the space 
charge limited regime is similar but the equations and 
initial conditions are slightly different. It is interesting 
to investigate whether the scaling relationship with ap- 

3/2 r 

plied voltage is different from the classical V g law of 
Child-Langmuir with the new boundary condition. The 
variation of the scaled critical current J c = J c / J s with 
U g = Vg/V s is shown in Fig. g] Clearly, J c - U% lo- 
cally with the exponent in the range 1.1 < a < 1.5 
and converging to 1.5 as U„ is increased. Note that 




FIG. 3. The critical perveance A c is plotted against E. The 
dashed lines are the classical prediction while the solid lines 
are obtained by numerically solving Eqns. (Q and ([5]). Inset 
shows the value of A c at small E for cj> g — 0.5. 



Thus, as V g and D are increased, J c obeys classical scal- 
ing. Note that the a = 1/2 behaviour reported in Ref. 
is not seen in our computations with the first order WKB 
boundary conditions. It is in fact possible to show^ that 
the plane wave boundary condition predicts a = 1/2. To 

see this, note that J ~ Ug^ 2 n(l) where n(l) = q 2 {l) and 



ip(x) = q(x)^/n^e l9 ( x \ Since 
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and q'(l) — for a plane wave boundary condition, 
dn/dU\jj = jj g = 0. Thus, n(l) is independent of U g . It 

follows that J ~ Vg 1 ^ 2 / D 4 only if the transmitted wave- 
function is assumed to be a plane wave. 




U g = Vg/V s - V g D 2 while J c 
J C D 5 ~ V n a D 2a or 

c 9 



Jc/Js - J C D 5 . Thus 
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FIG. 4. (Color online) The scaled critical current J c follows a 
local power law relationship with the scaled applied voltage: 
J c ~ Ug. Boxes mark the numerically obtained values while 
the short and long dashed lines are the fits at small and large 
U g respectively. At small U g , a ~ 1.12 while for large U g , 
a -> 1.5. 
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Note that the smallest valued of the scaled potential 
considered is U g = 2. At D = lOnm, this translates 
to an applied potential of 0.75mV while at D = 70nm, 
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it is 0.015mV. On the higher side where classical scal- 
ing holds, the maximum scaled potential considered is 
U g = 150000. At D = lOnm, the maximum V g = 57V, 
while at D = 70nm, the maximum applied voltage con- 
sidered is V g — 1.16V. Beyond these voltages, classical 
scaling should be applicable. These results are consistent 
with recent experiments on aluminium 15 and graphene^ 
nanogaps. In case of aluminium for instance (see fig. 3 
of Ref. [H), the 70nm gap shows good agreement with 
the Fowler-Nordheim emission law for V g J$ 6V. Above 
this (see Fig. 4 of Ref. 15), the current appears to be 
space charge limited with the exponent a ~ 1.5 which 
is consistent with the first order WKB results presented 
here. 

In case of graphenei^, the space charge limited regime 
occurs at smaller voltages possibly on account of Klein 
tunneling^ 7 -. For a few hundred nanometer gap, the ex- 
ponent a averaged over several gate voltages, takes the 
classical space charge limited value of 1.5 for V g gl 3V. 
This is again consistent with the results presented here. 

In conclusion, we have shown that a first order WKB 
approximation for the transmitted wave gives semiclas- 
sically consistent results which are significantly different 
from those obtained using a plane transmitted wave. The 
maximum transmitted current is found to be smaller than 
the classical current for most injection energies and may 
exceed only in the quantum limit of small E. Finally, at 
zero injection energy, the limiting current density follows 



a local scaling law J c ~ / D 5 ~ 2a with a — > 3/2 for 
large V g and D. Importantly, our result using first or- 
der WKB approximation^ underscores the possible pa- 
rameter regime where nonclassical space charge limited 
behaviour may be observed experimentally. 
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